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Abstract. The bead process is the particle system defined on parallel lines, with underlying 
measure giving constant weight to all configurations in which particles on neighbouring lines 
interlace, and zero weight otherwise. Motivated by the statistical mechanical model of the 
tiling of an abc-hexagon by three species of rhombi, a finitized version of the bead process 
is defined. The corresponding joint distribution can be realized as an eigenvalue probability 
density function for a sequence of random matrices. The finitized bead process is determinan- 
tal, and we give the correlation kernel in terms of Jacobi polynomials. Two scaling limits are 
considered: a global limit in which the spacing between lines goes to zero, and a certain bulk 
scaling limit. In the global limit the shape of the support of the particles is determined, while 
in the bulk scaling limit the bead process kernel of Boutillier is reclaimed, after approriate 
identification of the anisotropy parameter therein. 



1. Introduction 

Consider a particle system confined to equally spaced lines that are parallel to the y-axis and 
pass through the points x — k, k G Z+, and let this latter integer label the lines. For each line k 
place point particles uniformly at random, with the constraint that their positions interlace with 
the positions of the particles on lines k±l. It has been shown by Boutillier that such uniformly 
distributed interlaced configurations — referred to as a bead process — are naturally generated 
within dimer models. Furthermore, it is shown in 6J that the bead process is an example of a 
determinantal point process: the fc-point correlation function is equal to a fc x fc determinant 
with entries independent of k. If all the particles are on the same line, this correlation function is 
precisely that for the eigenvalues in the bulk of a random matrix ensemble with unitary symmetry, 
for example the Gaussian unitary ensemble of complex Hermitian matrices (GUE). 

If some of the particles are on different lines, and the anisotropy parameter (see below for this 
notion) is zero, the correlations between particles on different lines coincide with the correlation 
between bulk eigenvalues in the GUE minor process [9j. The latter is the multi-species particle 
system formed by the eigenvalues of a GUE matrix and its successive minors. Let Af„ denote an 
nx n GUE matrix, and form the principal minors Mt, t = l,...,rt— 1 as the txt GUE matrices 
corresponding to the top t x i block of Af„. Let the eigenvalues of Mt be denoted by {x'^p}j=i....,t 

with x^t^ < ■■■ < a;^*^ . The successive minors Mt+i, Mt have the fundamental property that 
their eigenvalues interlace, 

(1.1) < 4^ < 4^''' < ■ ■< 4^'^ < 4^ < 4^'^ 

for < = 1, . . . , n — 1. It turns out that for nx n GUE matrices [2], the joint distribution of all the 
^n{n + 1) eigenvalues ^s^i{x'f^}j=^,---,'^ successive minors is proportional to 



j=l l<j<k<n 
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where x{x^^\ x'"^'^^^) denotes the interlacing (|l.ip between neighbouring species. Regarding 
species (s) as occuring on line s, one sees that the eigenvalues to the left of line n all occur 
uniformly within interlaced regions, precisely as required by the bead process. 

The main objective of this paper is to construct and analyze a finitization of the bead process. 
Like in [B], we are motivated by a statistical mechanical model — in our case the tiling of an 
a5c-hexagon by three species of rhombi. However unlike in [B] this statistical mechanical model 
plays no essential role in the ensuing analysis. We will see in Section 2 that the rhombi tiling 
of an a5c-hexagon leads to a particle system defined on the segments < y < 1 of the lines at 
X = 1,2, . . . ,p + q — 1 (p < q) in the xy-plane. The number of particles, r{t) say, on line x = t 
(to be referred henceforth as line t) is given by 

( t, t<p 

(1.3) r{t) p, p<t<q 

[ p + q-t, q<t. 

Denote the coordinates on line t by a;*-*^ := {a;j*'}j=i_...^r(t)i ^-nd require that the particles be 
ordered Xr{t) < ■ ■ ■ < xi. With xix^'^\ a;*^*+^^) denoting the interlacing p.l|) between particles on 
neighbouring lines the particle system — which is our finitization of the bead process — is specified 
by the joint probability density function (PDF) on the line segments proportional to 

p— 1 q—l 

(1.4) nx(^^'^^^''''^)nx(^^^^^*+'^u{o}) n xix^'\x('+'^u{o,i}). 

t—l t—p t—q 

The particles in the tiling problem are restricted to lattice sites. To obtain our finitization 
of the bead process requires taking a continuum limit in which the lengths of the vertical side 
of the afoc-hexagon is taken to infinity. This is done in Section 3. There we also compute the 
marginal distribution of (|1.4p on any one line, and we use that knowledge, together with some 
random matrix theory, to compute the shape formed by the support of the particles in the limit 
that the lines form a continuum in the horizontal direction. 

We know that the PDF (|1.4p can be realized as a joint eigenvalue PDF as well as a distribution 
relating to both queues [2j and tilings [11]. Likewise, in Section 4 we show that the PDF (|1.4p . 
which can be obtained as the continuum limit of a tiling problem, can also be realized as an 
eigenvalue PDF for a certain sequence of random matrices. 

Section 5 is devoted to the computation of the general rt-point correlation function for n parti- 
cles on arbitrary lines in our finitized bead process. The joint PDF p.4p is a determinantal point 
process, meaning that the n-point correlation function can be specified as an n x n determinant 
with entries independent of n. To compute the entries, referred to as the correlation kernel, we 
use a method based on a formulation in terms of the diagonalization of certain matrices, due to 
Borodin and collaborators 3 . In Section 6 we take up the task of computing the bulk scaling of 
the correlation function. This involves us choosing the origin at the midpoint of a collection lines 
spaced 0(1) apart, scaling the distances on the lines so that the mean density in the ncigbour- 
hood of these points is unity, then taking the limit that the number of particles goes to infinity. 
We reclaim the bead process correlation kernel obtained by Boutillier [B], after identification 
of the anisotropy parameter therein with a parameter defined in terms of the continuum line 
number, and an anisotropy parameter relating to our underlying a5c-hexagon. 

2. Discrete model 

Construct a square grid from the lines x — ±j, y = ±fc, {j,k g Z) in the xy-plane. On this 
grid construct continuous lattice paths which consist of straight line segments connecting the 
lattice point {x,y) to {x + l,y ±1) for one of ±. Such a lattice path, starting at {x,y) and 
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finishing at {x + N,y+p) say, can be thought of as the space time trajectory of a random walker 
confined to the y-axis. The random walker starts at y, and finishes at y+p after N steps, each of 
which are one unit down, or one unit up. We are interested in n such lattice paths, conditioned 
never to intersect, starting at (0, 2i) and ending at (p + q, —p + q + 2i) for i = 0, . . . , n — 1. 
Such configurations are in bijection with domino tilings of a hexagon as seen in Figure [TJ It 
can be shown, [T7] that the positions of the lattice paths along the lines, i.e. the red particles in 
Figure [U form a determinantal process. But our interest is not in the particle process defined by 
the paths, which is further studied in [231 [H] but rather the process defined by the empty sites 
(holes) between the paths, which in Figure [Dare the blue particles. It suffices to assume p < q, 
as the other case reduces to this under refiection. 

Along line t there will be r{t) blue particles, where r{t) is given by (|1.3p . Let {xj*-'}^^!^ ,,(«)> 

with Xj.(t} < ■ ■ ■ < xi he their corresponding positions. The highest possible position for cc^*'' is 

(*) 

r(t) 



a{t) and the lowest possible position for is b{t) where 




With t in these equations regarded as a continuous parameter, < t < p + g, these lines form 
the four non- vertical sides of the hexagon in Figure [T] 

Let us introduce some fixed (virtual) particles by setting 

(2.2) - 2 = 6(t) - 2 (t = 0, . . . ,p - 1) 
and 

(2.3) xf ■.^2(n + q)-t^a{t)^1 {t = q, . . . ,p + q ~ 1). 

Then it is a consequence of the combinatorics of this model, as can easily be seen in Figure [T] 
that the blue particles fulfill certain interlacing requirements, 

(t) (t+i 



(2.4) xn,<x^Uxf^, 



for t = 0, . . . , q — 1 and 

(2.5) xn,<4^'Uxf\ 

for t = q, . . . , p + q — I. The probability of Uj^q {a;^*^}j=i_...^j.(t)} is therefore proportional to 

q—l p+q-l 

(2.6) n*(^^*^^^*'''^) n x(xW,x(*+i)). 

t=0 t=q 

Here x{x,y) refers to the ordering between sets x = {xi, . . . ,Xr} and y = {yi, . . . ,yr} of the 
same cardinality, 

Xr < yr < ■ ■ ■ < xi < yi. 
Note that ()2.6p is a discrete version of the joint PDF for our finitization of the bead process. 

The probability p.6p can be written in terms of determinants. Thus, using the identity 

(2.7) = det[Xy,>a:Ji,fc=l,...,r 

(see e.g. O lemma 1]), and writing (f>{x,y) := Xy>x, (|2.6p reads 

(2.8) ndet[</>(xf,xf^))]::^)x J] -f^''%%,^ 

t=0 t=q 
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Figure 1 . In the text the spacing between lattice paths is 2 units, while in the x- 
direction the spacing between lines is 1 unit, but in the figure these lengths have 
been scaled so that the six internal corner angles of the hexagon are equal. The 
three independent side lengths-left, bottom left, top left — are in correspondence 
with the number of walkers n, and the parameters p and q of the text respectively. 



One way to think about this expression is that the blue particles themselves are the positions 
of some dual random walks, with step probability (p. By the Lindstrom-Gessel-Viennot method, 
the measure on configurations is given exactly by (j2.8p . 

A natural question to ask then is what is the probability distribution of the blue particles 
along time t. The idea is to compute the number of possible configurations to the right and to 
the left of t and divide by the total number of possible tilings. For this we need the following 
lemma. 

Lemma 2.1. Let t < p. Given some configuration x^^\ the number of configurations to the left 
of line t, i.e. 

(2.9) ndet[0(xf\4^+^))],^+l, 
is 

(2.10) Ht{xi,...,xt) = ctl^{xi,...,xt). 
where A is the Vandermonde determinant and 

(2.11) ct ^ 



2t(t-i)/2nt-;fc, 
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Proof. Induction over t. For t = 1 it is indeed true that Hi{x) = 1. Assuming the statement is 
true for t, consider t+1. By summing the distinct variables afong distinct rows of the determinant 
we have 
(2.12) 

Ht+i{xi, . . . , xt+i) = 



Xi<yi<Xi+i 



1 2/1 



(2.13) 



1 yt 

Po{xi) -Pq{x2) 



yt 



Pl{xi) -Pl{x2) 



where pkix) = YTy=2Mi'^y)^ fo'" ^ ^ven and pk{x) = I]yl2M + 1)'' for x odd, M is an 
arbitrary big negative number. It can be checked that pk is a polynomial and that Pk{x) = 
x^^^ /2{k + 1) + 0{x^). Performing column operations in the determinant removes dependence 
on all but the highest order coefficient in pk and gives 



Po{xt) - PQ{xt+i) pi{xt) - pi{xt+i) 



Pt-l{xi) -Pt~l{x2) 



Pt-l{xt) -pt-l{xt+l) 



(2.14) 



Uxi-X2) \{xl-xl) 



■k;{x\ x\) 



\{xt-xt+i) \{x1-xl^^) 



Pull out constants and border the matrix to make it t + 1 by t + 1 to obtain 



(2.15) 



2H\ 



Xi — X2 



xt- Xt+1 xj - xf^^ 



x\ — 



t+1 



1 Xt+l 

Row operations now complete the proof. 



□ 



Theorem 2.2. The probability that at time t the blue particles are at positions xi > ■ ■ ■ > Xj.(t) 
is 

r{t) 



(2.16) 
where 
(2.17) 



Ptixi, 



r(t)) 



Zt ^^'^{xi,---,Xr(t))Y[Mxi) 



1=1 



\q-t\ \p-t\ 

Mx) - Yl ia{t) + 2k-x) Y[ix^ bit) + 2k) 

k=l fe=l 



and Zt is some normalizing constant. 

Note that this is called the Hahn ensemble in the literature, see for example [17] . 

Proof. Case t < p: The area to the left of t can be tiled in Ht{xi, . . . , xt) ways and the area to 
the right of t can be tiled in Hp+q^t{x^q+t+i, ■ ■ ■ ,xi, . . . ,Xt, ■ ■ ■ , Xp) ways where we introduce 
virtual particles X-q+t+i = a{t) + 2{q — i — i + 1), for z = 1, . . . , g — i and Xt+t = b(t) — 2i for 
z = 1, . . . , p — t (see Figured]) Then 



(2.18) ptixi 



r(t) 



) = const ^Ht{xi, . . . ,Xt)Hp+q^t{X-q+t+l, ..■,Xl,...,Xt,...,Xp) 
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Figure 2. Construction of virtual particles as used in the case t < p oi the 
proof of Theorem 12.21 The blue particles can be seen as the positions of some 
dual random walkers, leading to a probability distribution proportional to the 
expression in ()2.8p . 



which proves the theorem in this case since the part of the Vandermonde determinant that has 
to do with the virtual particles is exactly ft- 

Case p < t < q: The area to the left of t can be tiled in Ht{xi, . . . , Xp, . . . , Xt) ways where 
we introduce virtual particles — b(t) — 2i for i — 1, . . . , t — p. The area to the right 

of t can be tiled in exactly Hp^q^t{x-q+t+ii ■ ■ ■ ,xi, . . . ,Xp) ways where the virtual particles 
x^q^t+i — + 2i for i — I, . . . , q — t. The theorem then follows just like in the first case. 

Case q < t. The area to the left of t can be tiled in Ht{xq-t+i, . . . , xi, . . . , Xp) ways with 
virtual particles Xq^t+i = ait) + q~t~2i for i = 1, . . . , t — q and Xp+q-t+i = b(t) — 2i for i = 1, 
. . . , t — p. Again the theorem follows. □ 

In [in] it is shown that the positions of all the blue particles — with the probability measure 
defined by (|2.6p — is a determinantal point process and its kernel is computed in terms of the 
Hahn polynomials. An alternative proof of Proposition lFTl would be to take their expression and 
applying the known limit formula 

(2.19) hm Q,,{Nx, a, /3, N) = "J.^^ ' 

N^oo P,?'' (1) 

where Qn and P"'^ is the n:th Hahn polynomial respectively Jacobi polynomial. Our motivation 
for proving Proposition 15. II from first principles in section [5] is that our proof is more direct and 
intuitive. Also, the idea of virtual particles in fixed positions, as illustrated in the previous proof, 
has not been previously used to our knowledge. That idea also gives a simple proof, see [35], of 
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a result conjectured in [T^ and proved in |23j, that the GUE minor process can be obtained as 
a scaling limit of lozenge tilings of a hexagon close to the boundary. 

3. Continuous model 

We are interested in the measure corresponding to the probability (|2.8p under the rescaling 
Xj*^ 1-^ Xj'^ jn when n ^ oo while p and q are kept fixed. It is easy to see that this is exactly the 
measure corresponding to the PDF (|1.4p . We have thus derived our finitized bead process as the 
continuum limit of a tiling model. In fact a semi-continuous lattice paths model (see Figure |3j 
equivalent to our finitized bead process was proposed in [1 to model the bus transportation 
system in Cuernavaca (Mexico). In that context, for p < t < q the particle configuration 
,,(*) is the time at which bus number j arrives at bus stop labelled t. 

For future reference, we note that analogous to the equality between (I2.6P and (|2.8p . use 
of ()2.7p allows (11. 4p to be rewritten as 

q-l p+g— 1 

(3.1) ndet[0(a.f\xr^^)]-?^^^x J] det[0(.f), ^^^^ )]:;£„ 



which apart from the absence of the t — term in the first product is formally identical to 
The variables take values on different sets though, with those in (|2.8p being confined to lattice 
points, and furthermore in (j3.ip the virtual particles are now specified by 

(3.2) xft+D^O (i= l,...,p-l) 4*^ = 1 {t = q,...,p + q-l) 



instead of UM- 

Rescaling the measure in Theorem 12.21 gives the following (cf. [H eq. (10)]). 

Proposition 3.1. The projection of the PDF |i.^[ ) to x*^*) is for xi > ■ ■ ■ > Xr(t) the density 

r(t) 

(3.3) Pt{xi,...,Xr{t)) = Z^'^ ^'^{Xl, ■ . ■ ,Xr{t))W_ft{Xi) 

1=1 

where 

(3.4) /t(x) = (l-x)l«-*lxlf-*l 
and Zt is some normalizing constant. 

The PDF (|3.3p is familiar in random matrix theory as the eigenvalue probability density 
function for the Jacobi unitary ensemble (see e.g. [8]). Due to this interpretation, it has appeared 
in a number of previous studies, and some properties of relevance to the bead process are known. 
Two of these are the support of the limiting density, and the limiting density itself. To specify 
these results, for notational convenience let us write r{t) = n and 

(3.5) ft{x)^x''{l-xf, 

and consider p.3p in the case that a — an, [3 — bn and n — > oo. Then we know from [25|, [71 [H] 
that the limiting support is the interval [c, d], where c and d are specified in terms of a, h by 

(2 + a + 6)2 + (a2„^2) 
c + d = - 



(2 + a + 6)2 

(3.6) d-c^ ^/{c + df -Acd, 

and in this interval 



(3.7) P(i)(y) = hm -p(i)(y) 



^ 2 + a + b^/{y~c){d-y) 

%Z 2^ y(i-y) ■ 
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Figure 3. In this lattice paths model introduced in H], up/right comers corre- 
spond to arrival times at stations 2,3,4,5. Regarding all the up/right corners 
as particles our finitized bead process results. 



To make use of these results, let us introduce a scaled label S for line t, and a scaled measure 
of the difference in the left bottom and left top side measures of the hexagon by 

(3.8) S^t/p, k^{q-p)/p 

respectively. We can then analyze the limit of the finitized bead process in which t,p,q — > oo 
with S, k fixed. Note that the scaled line label can take on values 

(3.9) 0<S'<2 + fc. 

Straightforward application of (|3.6p and (|3.7p reveals the following limit theorem for the finitized 
bead process. 

Theorem 3.2. In the above specified setting, the support of the head process along scaled line S 
is [65, cs] with 



Sk 1 2^5'(fc + l)(fc + 2-S') 

~ (fc + 2)2 ^ 1+2 (fc + 2)2 

Sk 1 2^S{k + l){k + 2~S) 

^'■''^ ^^^(fcT^ + fcT^ + (kTW ■ 

The normalized density is given by Jg. 7[ ) with c = cs, d — ds and (a, b) given by 
(3.11) {^^, o ' {S-l,k+l~S) 



S ' S J' ' ' " \2 + k- S' 2 + k- S 

for 0<5<1, 1<S'<1 + ^ o,nd l + k<S<2 + k respectively. 

We note that at S" = 0, cs = ds = l/(fc + 2), and similarly a,t S = k + 2, cs = ds = 
{k + l)/(fc + 2). The degeneration of the support to a single point at the left and right sides is 
in keeping with the number of particles in the bead process starting off at one at these sides. 
In p.4p . when t = p the exponent a in p.Sp vanishes. Thinking of the factor x" as a repulsion 
from X — 0, with this term not present we would expect that the left support to be its smallest 
allowed value, cs = 0, which is indeed the case in (|3.10p with 5 = 1. Similarly, in (|3.4p with 
t ^ q the exponent /3 in (|3.5p vanishes, and correspondingly in p.lOp ^5 = 1 for S* = 1 + fc. 
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A particular realization of the boundary of support as implied by (|3.10p is plotted in Figure 

m 

4. Realization as a random matrix PDF 

We have seen that the joint PDF (|1.4|1 defining our finitization of the bead process arises 
naturally in the context of the tiling of an a6c-hexagon, or equivalently by the consideration of 
non-intersecting paths. Here it will be shown that (|1.4p can be obtained as the joint eigenvalue 
PDF of a sequence of random matrices. 

Our construction is based on theory related to certain random corank 1 projections contained 
in [21 US], which will now be revised. Let 

M ^ nAU, U = I-xx'^ 

where 

A = diag((al)^^(a2)^^...,(a„)^"). 

Here the notation (a)^' means a is repeated p times, and it is assumed ai > 02 > ■ • ■ > a„, while 
if is a normalized complex Gaussian vector of the same number of rows as A. The eigenvalues Ui 
of A occur in M with multiplicity Si — 1. Zero is also an eigenvalue of M. The remaining n — 1 
eigenvalues of M occur at the zeros of the random rational function 

n 

where (gi, . . . , (?„) has the Dirichlet distribution D[si, . . . , s„]. With these n — I eigenvalues 
denoted Ai > • • • > A„_i, it follows from this latter fact that their joint distribution is equal to 
(4.2) 



r(si H h Sn) Ill<j<k<n-li^J ^ ^k) 



' n — 1 n 



j ~ '^pl'^'' ^ I Xai>Ai>a2> - >A„_i>a„ 



j=l p=l 



(as done in (|1.2p . the indicator function will be abbreviated x(A,a)). 

After this revision, we begin the construction by forming Mi = H^iH, where Ai = diag((0)^, (1)'). 
It follows from the above that Mi has one eigenvalue Aj^^ different from and 1 satisfying 

< Ai^^ < 1, and this eigenvalue has PDF proportional to 

Now, for r = 2, . . . ,p inductively generate {x'f^}i=i^...,r as the eigenvalues different from and 

1 of the matrix 

Mr = UArU 

where 

Ar = diag((0)^-'-+\ Ar^\ . . . , xt-'\ (1)'-^+^) . 

It follows from (|4.2p that the PDF of {A^'^^}j=i,..._r, for given {A^'^^^^}j=i,..._r-i, is proportional 
to 

(4 3) x(a('^-i) -< ni<j (Ar - nLi(i - >^k^)''~''i>'k^T~'" 
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Forming the product of (|4.3p over r = 1, . . . ,p gives the joint PDF of U^^ilx'f^} which is therefore 
proportional to 

(4.4) n x(A('-i), AW) f[{X^^ - ) - Xi^^r^^. 

r—2 i<j k—1 

Next, for r — 1, ... ,17 — p, inductively generate {x[^^^'^}i=i^,,,^p as the eigenvalues different 
from and 1 of 

where 

Ap+r = diag(A(P+'--^\ . . . , A(^+-i), iiy-P-r+^) . 

We have from g^]) that the PDF of {X'f^'''^}j=i^,,,^p, for given {x'f^''^'^'>}j=i^,,,^p, is proportional 
to 

(4.5) ^(A(p+''-i),A(^+'')u{0}) ^^-</V' <Ju • 

The joint PDF of U^^]^{Aj*''} is obtained by multiplying the product of (|4.5p over r — 1, . . . , q~p 
by (|4.4p . It is therefore proportional to 

(4.6) nx(A('-i),AW) n x(A('-i\AWu{0})n(A(')-A(')). 

r—2 r— p+1 z<j 

The final step is to inductively generate {X'f^^''}i^i^,,,^p-r {r — 1, ... ,p — 1) as the eigenvalue 
different from of 

where 

A,+,=diag(Al''+'^-^\...,A(''+'-i)). 

According to g^, the PDF of {A^-''+'''}j=i,...,p_r for given ,^p-r+i is proportional 

to 

' ' ' '^nr<7-^(Ar-^^-A5^--))ncr^(i-Ar-)).--i- 

Forming the product over r — 1, . . . ,p — 1 and multiplying by the conditional PDF (j4.6p gives 
that the joint PDF of U^tV^ i^'j"^} is proportional to 

(4.8) n^(A^'"'^A^''^) n x(A('^-^),AMu{0}) n x(A(''-i\a('-)u{0,1}) 

r—2 r— p+1 r— g+1 



and thus to the finite bead process ()1.4p . 

Graphical output from the implementation of this construction for particular p and q is given 
in Figure m 
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Figure 4. A configuration of the finite bead process with p = 4, q = 12, 
generated from the zeros of a sequence of random rational functions (|4.H) , cor- 
responding to the eigenvalues of the sequence of random matrices described in 
the text. The bounding curves are the limiting shape as implied by (|3.10p . 



5. Correlation functions 



It was shown in [T^ that the GUE minor process as specified by the joint PDF (|1.2p is 
determinantal, with the correlations between n eigenvalues chosen from any of the minors being 
given by an n X n determinant. The functional form of the entries of the determinant is of 
particular interest. Thus we have from [19j that 



(5.1) 
with 



P{n){{{Sj,yj)}j 



(5.2) 



K 



GUEm 



iis,x), {t,y)) 



Hs-k{x)Ht-k{y) 



k=l 



-)t-k 



k— — oc 



s > t 



s > t 



(here Hj{x) denotes the Hermite polynomials of degree j). We seek an analogous formula for 
the correlations of our finitized bead process. For this we make us of a linear algebra method 
developed by Borodin and collaborators [H [3] , although an alternative strategy would be the 
method of Nagao and Forrester [2j . Either method can be used to calculate (|5.2p , and various 
generalizations [^[TT]. 

Formalism. We take line j, which is the segment in the a;y-plane {j,t), < t < 1, and replace 
it with a discretization Mj. In the limit that \A4j\ — > cx) it is required that the discretization 



be dense on the segment. On the sets A4j U {0} {j — 1, . 
{a;p^}i=i,...j+i with x'fh = 0; on the sets Mp-i+j {j = 1, . 



. ,p — 1) distribute j + 1 points 
. ,q — p + 1) distribute p points 
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[x. 



(P-I+J) 



}^ 



{a:i^+^'^},=o,...,p-, with 4'+^") = 1. 



and on the sets Mq+j U {1} (j — 1, ... ,p — 1) distribute p + 1 ~ j points 



On the configuration of points X 



} define a (possibly signed) measure 



(5.3) 



1 

-ndet[T4^,( 



C 



t=l 



p+q~2 

n det[M^z(, 

t=q 



Jt) „(t+l)^lKt) 
' i-i ''ij=i 



for some functions {T^/^ : A^i x Mi+i — > R}i=i,...^p+g_i. With ^^((y, a;) = (j){y,x) = Xx>y, this 
corresponds to a discretization of the PDF (|3.ip for our finitized bead process. 

Next introduce the \Mi\ x |A^i+i| matrices Wi = \Wi{xi,yj)\xi£MuyjeMi+i- Introduce too 
the p X |A^;+i| matrices Ei (/ = 0, . . . ,p — 1) with entries in row / + 1 aU Is and aU other entries 
0, and the Mp^q-i x p matrices Fi {I ~ 1, . . . ,p) with entries in column / all Is and all other 



entries 0. In terms of these matrices, with A4 
\M\ X \M\ matrix L by 



{1, . . . ,p} U Ml U • • • U Mp+q-i, define the 






Eo 


El 


E2 


Ep-1 














-Wi 























-W2 ■■• 



















































-Wp ■■. 





Fp-1 
F2 

L i^i 





-Wq+l 





-Wp+q-2 





for block matrices of appropriate size. Note that the rows and columns of L are labelled by 
the set M. 

For Y C M, introduce the notation Ly to denote the restriction of L to the corresponding 
rows and columns. Then, in accordance with the general theory of L-ensembles [4^, our matrix 
L has been constructed so that the measure ()5.3p is equal to 



(5.4) 



det L^i p^tjx 
det(l* + L) 



where 1* is the \A4\ x \A4\ identity matrix with the first p ones set to zero. The significance of 
this structure is the general fact that the correlation function for particles a,tYcM/{l,...,p} 
is then given by 



(5.5) 



p(y) =deti^y, 



K = 1* - (1* +L)-i 



M/{l,...,p}- 



Our aim th.cn is to evaluate K in the continuum limit, with W^^(y,x) — Xy<x 

independant of 

i, and to do this we must turn (1* + L)^^ into a form that can be computed explicitly. As a 
start, write L in the structured form 



(5.6) 



L 




C 



B 
D - 1 
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where 

(5.7) B = [En, El,..., Ep^i, 0,0,..., 0] 

(5.8) C = [0, 0, . . . , 0, Fp_i, Fp_2, . . . , Fif 
We know that g] 

(5.9) K = 1* - D-^ + D-^CM-^BD-'^ 
with AI :— BD^^C and furthermore 

(5.10) D-' = l + [l^[,,,)].,,=i,...,p+,-i 
for 



(5.11) W< 



W,...Wj.i, Kj 
0, I > J. 



>From (|5.10p we compute that the m-th member of the block row vector BD ^ is equal to 

m — 1 

(5.12) E^.i + Ek-iW[k,m) 

k=l 

for 1 < m < p and 

p 

(5.13) ^£;fc_iTV[fc,™) 

for p+l<m<p + q— 1. Similarly the m-th member of the block column vector D^^C is equal 
to 

p-i 

(5.14) J2^l--^i+k)Pp-k 
for I < m < q and 

p+g— 1— m 

(5.15) Fp+q-m + ^ W[jyi^jn+k)Fp+q-m-k 

fc=l 

for q + l<m<p + q— 1. 

Explicit formulas. Our aim is to find square, p x p matrices Bq, Cq such that: BD~^ — Bq^ 
for some block row vector $ — [$^, . . . , D~^C = ^Cq for some block column vector 

— [^^, . . . , and BqCq = M. Finding these square matrices will allow us to compute 

K explicitly, since if the above equalities hold, then 

(5.16) K = 1* - D-^ + •i'<i> 
We introduce the notation 

(5.17) {a*b){y,x)= / a{y , z)b{z , x)dz 

Jo 

and note that, in the continuum limit Mi — > [0, 1], (i — I, . . . ,p + q — 1), and with the lattice 
points weighted by their mean spacing, (|5.1ip becomes 



(5.18) W[,,,){y,x) 



(Wi * • • • * Wj_i)(y,a;), i<j 
0, i>j 
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Independent of i, Wi{y, x) = Xv<x and so 

(5.19) Wi,^j){y,x) = J—J^7Jy Xy<x{x ~ yy-'-^ 

With the entry of K corresponding to hne s, position y for the row, and hne t, position x for the 
column, denoted K{s, y; t, x) it foUows from (|5.10p , (|5.16p and (|5.19p that 

1 ^ 

(5.20) K{s, y- 1, x) = - _ (x - yf-'-'xyKx + ^(*^),,i(<i>*)z,. 

^ 1=1 

so our task is to calculate the quantities {'4)^)y^i, {^*)i^x- 

We begin by looking at the m-th member of the block row vector BD~^ for m < p as shown 
in ([51^ . Using ((5715)) we see that 

Cm~l \ m—l 

where the convention ^ = for a e Z<o is to be used for m < I. Next define the family of 
polynomials 

(5.22) Pi''^''\x) = Pt''\l - 2x) = { n\x-{l-x)>>dx^^ ^ ' ^ - ^ 

0, n < 



These polynomials have orthogonality 

(5.23) x'^il - x)'Pj''-''\x)Pt-'\x)dx = AAj'^''^<5,fc 



where 

1 (n + a)!(n + 6)! 

2n + a + 6+l n!(n + a + 6)! 



(5.24) A/'i'^''') 



and as such, as indicated in (|5.22p . are the Jacobi polynomials P^'^'^ modified to be orthogonal 
on (0,1). 

Expressing (|5.2ip as a sum of these polynomials for certain a=p — m, b = q — m gives 



X 



m-l 1 ™ p(p-^,'l-"^)( -) .1 

y ^n-j K) / ^m-l ^p-™.^ _ t)9-™p(f-;"'«-™)(t)di 



(m-l)l (m-l)l'^ i^(p-m,q^m) J 

j=l •'^m-j 

where the second equality follows using ()5.22p and integrating by parts m ^ j times. If we let 

A/p_„ '^'ip-ny.{n-iy. Jo 

and 

/ -I'.p+m (P + g ~ ~ p(p-m,g-m), X 

(5.27) ($'")„,,= <( ^ (g_n)! - 

) n > m 
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then 
(5.28) 

and therefore 
(5.29) 



{m-iy. 



k=l 



Consider next the m = {p + j)-th member of the block row vector BD ^ as shown in (j5.13p . 
Using (|5.19p we see that 



(5.30) E^^-i%-.p+^) 
To start with, we look at the m = p case of (|5.28p 
(5.31) 



p+j-i 



{p+j-iy. 



ip-iy. 



E(^o).,„pi°r^^(-) 



We introduce the operator J[f{x)] — f{t)dt and note that 



(5.32) 



Applying to both sides of (|5.3ip gives 



(5.33) 

If we let 
(5.34) 



X' 



{p+j-iy. ^i(j-i)!yo 



{x-ty-'f{t)dt 



(c&P+J"), 



then it follows from ([QUI and (fOi]) that 



(5.35) 



J2Ek-iWik,p+,)=B„^P+^ 



k=l 



and this coupled with l|5.29p gives us our desired BD ^ = Bq^. 

Now that we have defined Bq, we can set about finding what C'o must be, using the equation 

(5.36) BoC'o = M 

First, using f^T]) . and (jSH})) we have 



(5.37) 



M = 



j=i \k=i ) 



The formula (j5.30p for the inner sum implies 
(5.38) (M)j,fe 



1 



(p + g + 1 - j - fc) 
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To proceed further requires the first of Jacobi polynomial identities 
d 



(5.39) ^ (^''■^"'''''^(^^) " {n + a)x''-^P^''-^'''+^\x) 

(5.40) ^((i_^)^pK'')(a;)) = ~{n + b)il-xf-'Pi^+'^''-'Hx), 

which can be verified from (|5.22p (the second will be of use subsequently). Setting j — q + I — 
m,l = n and a; = 1 in (|5.33p gives 

(5-41) r^—T T=E(So)„,. f'^^7^^Pi":r\t)dt 

[p + q+l-n-m) ^ (q - my. p " 

But it follows from ((5351l that 

(5.42) p(°-'^)(:r) = ^^^[x^ pi^^-^) [x)), a e Z>o, 

[n + ay. 

so the integral in ()5.4ip can be reduced by integration by parts, giving 

(5 43) 1 ^ (^o)«,fc f' t^-H^-tr-"' ^^ 

(p + q + l-n-m) (^ - «)! Jo {q - k)l 

Recalling (|5.36p and ()5.38p we thus have 

We now look at the m = (q + j)-th member of the block column vector D~^C as shown in 
((57T^ . We see that 

(5.45) ^ -> 



Expressing this as a sum of the polynomials from (15.22^ for certain a = j + q — p,b = j gives 



{p-j-ny {p-j-ny ^ ^f^^+1-p^3) 



(5.46) = ^(Co) 



(1 - t)p-j'"" tJ'+«-?'(i - ty p^^/_i^-'''' {t)dt 
(q-ky piitT'\y) 



We remark that this equation can alternatively be derived from (|5.25l) by writing x = 1 ~ y, 
m— p+l^j,l = n, and noting the symmetry of the Jacobi polynomials Pn°''''\l — x) = 
{-l)"Pi'''''\x). It follows from (lOSll and (fO^ that by setting 



(5.47) {<!l'*'),.k = < ip-j-k)! M"*n"' 

k > p — j 

we obtain 

(5.48) Fp^,+ J2 W[,+,.,+,+k)Fp-j-,^^'^+^Co 
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Finally, we look at the m-th block of the column vector D as shown in (I5.14p . We see 
that 



(5.49) 



\k=l 



(1 _ yy+q-m-n 

{p + q — m — n)l 



With m = q this corresponds to j = case in (|5.46p , so we have 



(5.50) 



{p-n)l 



k=l 



After making a substitution 1 — y ^ w and applying , j ~ q — m, we end up with 



(5.51) 

Setting 
(5.52) 

gives us 

(5.53) 



(1 _ yY+q-m~n 

{p + q — m — n)\ 



(Co)i,„ {q~l)\ f' (^ - 2/)''-i-™ 5 



)v.i 



E \-'UJi,,n \-i -r j \~ a I p 



(z)(iz 



and this, along with (|5.48p gives us our desired D^^C ~ ^'Cq. 

The quantities i^'')y,i and ($*);, ^ in (lOOl) have now been specified in ([OT]) . ([F^i]) . (fCTf|) 
and ()5.52p , allowing us (after a little bit more calculation) to write down the explicit form of the 
entries for the correlation function in terms of Jacobi polynomials. 

Proposition 5.1. Let {t, x) refer to the position x on line t of the bead process. For the config- 
uration Y = Ur=i{(^»' ■^«)} correlation function is given by 



(5.54) 

where 

(5.55) 

Here 
(5.56) 



p{Y) ^ det[K{t„Xi;t.j,Xj)] 



K{s,y;t,x) = < 



asivMx) J2 Mt) 
=1 ^1 



s{y)h{x) J2 



asiy) 



l=-oc 



(_y)P-(l_y)9- 

il-y)i-^ 

1 



(t) 



s > t 



s < t 



< s <p 
p < s < q 
q < s < p + q 



(5.57) 



btix) 



(-1) 



p-t 



0<t<p 
p <t < q 
a;*-P(l - x)'-« q<t<p + q-l 
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(5.58) 



C 



(t) 



{p-l)l 

{p + q-t-iy. 
{p-iy. 



(5.59) 



(5.60) 



<t<p 
p <t<q 
q<t<p+q-l 

<t<p 
p <t <q 



Pil7-tl\-) q<t<P + q-l 



p-l 



0<t<p 
p <t <q 
q <t <p + q-l 



and a{s, t) = m\n[p + q — s,t^p\. 



Proof. The s>t case follows directly from (|5.20p by imputing the values given in (|5.27p , ()5.34p , 
((CT7)) and (|53^ . and using ([O^ . ((05)) and (|5^ . For s < t some further calculation is 
required, as in this case the first term in (|5.20p is non-zero. In this case we express this term as 
a sum of polynomials in x 



(5.61) 



^ OO 



(or similarly in y), with the aim to obtain the form in (|5.55p . It turns out that this expansion is 
such that 



(5.62) 



bt{x) ^ aiQl^'{x) = as{v)ht{x) ^ 



' C^^^ Q\'\y)Qf\x) 



1 — — 00 



At) jr{t) 



1 = 1 



(5.63) 6,(x)^a,gp>(x) = Y.i'^ly.i^ 

1=1 

For example, for p < s < t < q 

(5.64) E(*%.i(<i>*)' 



l^{s-i)\p^^ir''\x)pi;_-r~'\y) 



1=1 

We expand in y 
(5.65) 



1=1 



{t-iy. 



,r(s-p,q-s) 



(t-s-l) 
Using the orthogonality (|5.23|) 



- OO 

—TY,{x yY-'-'xy<. = (1 - y)"-' E "'^/^ 



(s-p,q-s) 



(y) 



(5.66) 7V;- ■ - ■ai = 



■(s—p,q — s) f (-^ 

(i-s-l)! 



i=0 



Using (|5.39p and integration by parts {t — s) times gives 
(5.67) 
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(5.68) , ' ,M -Vr'-'Xy<.-x'-'{l-yY~' y ^T^^ — — - 



1 — — 00 I 

(5.69) =J2i^%,ii<^%,,+asiyMx) ^ 

1 = 1 i=-oo ■^l 

as required. □ 



6. Bulk scaling 

In Section 3, a scaled limit of the bead process was considered in which the lines themselves 
formed a continuum. We computed the support of the density, and the density profile. In this 
section we focus attention on a different scaled limit — referred to as bulk scaling — in which the 
particles are on lines a finite distance (as measured by the line number difference) apart, with the 
interparticle spacing on the lines scaled to be unity. We require too that the particles be away 
from the boundary of support (thus the use of the term bulk scaling), which we do by requiring 
them to be in the neighbourhood of the midpoint of the support, {cs + ds)/2 in the notation of 
Theorem O 

Explicitly, the location t of the lines will be measured from a reference line of continuum label 
S (recall p.Sp . and thus t i-^ pS + t. With Xs {cs + ds)/'2 denoting the midpoint of the 
support on line 5, and 

(6.1) Us := P{i){Xs) 

where p(^i^ is the normalized density p.7p on line S the location x of the particles are to be chosen 
so that X = X + Xi/ {r{pS)us)- (The quantity r{t) denotes the number of particles on line t as 
specified by (|1.3p ). Our objective in this section is to compute the scaled correlation kernel 

(6.2) K{so,Y;to,X) := lim — i- — K (pS + so,y;pS + to, x) 



with 



(6.3) x^Xs + —^—, y = Xs + -^—, 

r{pS)us r{pS)us 

and thus according to (|5.54p the scaled correlation function 



(6.4) p{h,Xi;...;tr^,Xn) 
f 1 \ " 

= 1™ \ ( p{pS + h,Xs + Xi/{r{pS)us)]. . .\pS + tn,Xs + X^/{r{pS)us)). 

The strategy to be adopted is to make use of a known asymptotic expansion for the large n 
form of p^+a"./3+f'"(a;) with x such that the leading behaviour is oscillatory (outside this interval 
there is exponential decay) [20] . Based on experience with similar calculations [TOl [9l [11] we 
expect that this will show that to leading order the sums are Riemann sums, and so turn into 
integrals in the limit p ^ oo. 
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Proposition 6.1. Let a,b,x be given, and define the parameters A,p, 0,7 according to 









A 


= [a{x + 


l) + 6(a;- 


1)]^ 


-4(a + 6+l 




2e 


ip 




a{x + 1) 


+ 6(x- 1) 


+ iv 






+ a + 


b){l 




(1- 


f a + 6)(l- 








2(a4- 


1) 




(a + 6 + 


2)x - (3a - 


f 6 + 


2)-i^/^ 


(1- 


x){l 


f a4 






2(x- 1)(1 


+ a- 






2(6 + 


1) 




(a + 6 + 


2)x+ (a + 


36 + 


2) -ix/^ 


(1 + 


x){l- 


f a + 






2(x+ 1)(1 


+ a- 


^b) 



(6.5) 



For A < we /lawe the large n asymptotic expansion 



-TT < p < TT, 



-TT < < TT 



-TT < 7 < TT 



(6.6) P,("+""'^+''")(j;) 

2(6 



2(a+ 1) 



1) 



(l + 2:)(l + a + 6) 



Kf>+i)+f + 



(l-2;)(l + a + 6) 



KQ+l) + f + ^ 



(l-a;2)(a + 6 + l) 



X cos ( [n{a + 1) + a + + [n{b + I) + [3 + - {n + ^)p + 



1 



wa/id for general a, /? G R /or a, 6 > 0. As noted in [7j, results from [T31[S] imply that the 0{l/n) 
term holds uniformly in the parameters. 



Actually (|6.6p differs from the form reported in [SU], with our \/— A in the denominator of 
the first term on the RHS, whereas it is in the numerator of the corresponding term in |.20 , and 
furthermore some signs and factors of 2 in the cosine are in disagreement. One check is to exhibit 
the symmetry of the Jacobi polynomials 

(6.7) P^'^)(-x) = (-l)"p('^'^)(x) 

For this we examine the effect on the parameters (|6.5|) under the mappings 

(6.8) X I— > —X, ai-^ b, b^ a 

We see that A is unchanged, while 

(6.9) p i—t t: — p, I— > —7, 7 I— )■ ~9. 

Making the substitutions (16. 8p . (|6.9p in (|6.6p . along with a 1— > /3, /3 i— > a we see that indeed the 
RHS is consistent with (|6.7p . Another check is to specialize to the case = 6 = 0. Then with 
X = COS0, < (/) < TT we can check from (|6.5p that V— A — 2sni(j), p = 7r/2, 61 = 7r/2 — (^/2, 
7 = — '/>/2 and so 

1 \i/2 1 



(a,/3) 



(C0S( 



(-) 



■ COS 



((71 + (a + /? + 1)/2)6I - (a + l/2)7r/^ 



(sin(/./2)"+i/2(cos0/2)'3+i/2 

which agrees with the result in Szegos book |24j . We give our working in the Appendix. 

According to (|5.55p and (|5.59p the particular Jacobi polynomials appearing in the summation 
specifying K in (j6.2p depends on the range of values of the continuum line label S. Let us 
suppose that \ < S < k + 1, and so the number of particles on each line is p. In this case, and 
with So > to 

(6.10) K{pS + so.y;pS + to,x) = ^ c„(a;,y)Pi''"+*'"''"-*'')(x)Pi°"+^-''"-^°)(y) 

n=0 
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where 

(6.11) c„(x,y) = (1 - ,)^n-.o^.n+.o (2" + «^ + + l)(n + an + bn)\n\ 
and with w :— n/p^ 

(6.12) a:^'-l, 5=i±A^. 

w w 

The expression in the case sq < to is the same except that the range of the sum is now over 1 to 
oo instead of to p — 1. With x and y given by ()6.3p and r{pS) therein equal to p (recah ()1.3|) ). 
we want to replace the summand by its large n asymptotic form. Because of the form of x and 
y in (|6.3p . the parameters A,/9, 0, 7 in Proposition 16. 11 as they apply to the Jacobi polynomials 
in (|5.59p . all have expansions in inverse powers of 1/p. 

Lemma 6.2. Write 

X 

(6.13) X y-^ 1 ~ 2x with x — Xs H . 

pus 



The quantities h6.5\l have the large p expansion 

A = Ao + AiX/p + 0(l//) 
p^Po+PiX/p+0{l/p^) 

e = eo + eiXip + o{\/p') 

(6.14) 7 = 7o+7iX/p+0(l/p') 



where 



Ao =4 (a2(l - Xsf - Xs{2a{2 + 6)(1 - Xs) + b^Xs - 46(1 - Xs) ~ 4(1 - Xs))) 
Ai =— ((4 + 9)Xs -2- a^{\ - Xs) - a(2 + h){\ - 2Xs) - 26(1 - 2Xs)) 

Us 

g.po _ 2a(l -Xs)-2bX5 + iV^ 
4^(l + a + 6)Xs(l-X5) 
1 a - (a - 6)^5 



Pi = 



V^MsXs(l-^s) 

2a + 2(q + & + 2)Xs + ^V^^ 
4v/(a + l)(l + a + 6)X5 
1 a(l - X5) - (2 + b)Xs 



\/-Ao us-'^s 

4^(6+l)(l + a + 6)(l-X5) 



2a + 46 + 4 - 2(a + 5 + 2)Xs - «V^A(^ 

1 (2 + 6)Xs-a(l-X5)-2 
6.15 71 ^ -j TT^^ . 

V-Ao Ms(l-As) 

Proof. According to (|6.5p . the expansions of p,9,'^ depend on the expansion A, so the latter 
must be done first. For this, all that is required is to substitute ()6.13p in its definition and 
perform elementary manipulations. Substituting the expansion of A into the definitions of p, 9, 7 
the values of po,0o, 70 can be determined immediately. Making use of these, the stated values of 
pi , 6*1 , 71 then follow. □ 

We can now proceed to compute the asymptotic form of the summand in (j6.10p . 
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Proposition 6.3. Let x,y be as in i6.3\) . and let wq < w < 1 correspond to A < 0. For large p, 
and to leading order in contribution to the summation, we can replace 

c„(x,y)Pi™+*'''''"-*«)(x)P^''"+'*'''''"-''«)(2/) 

by 



(6.16) '±±±±^n^o-..f-_ 



X Re 



y 

exp 



1 -X 



(X - Y)^/^>\ ( 2a(l - Xs) + 2bXs + ^^/^^^ 



so-*o' 



AusXs{l -Xs) J \ AXs{l - Xs) 

For Q < w < Wq, corresponding to A > 0, to the same order the summand can be replaced by 



Proof. Use of Stirling's formula gives 
(6.17) c„(x, y) - (1 - y)''"-*°x™+*°u;p(a + & + 2) 



{n + bn)'" (a + &+l)"(°+^+i)+5 



while we read off from (|6.6p that, provided A < and thus wq < w <1, 



(a + 1) 



Trn \^^A(a:)A(y) / 

„(a+l)+<±!l+M+i 



(1 + a + 6) 



(l + a + 6) 



-,„(b+i)_<io+i!i) + i 



where 



X (a + 6 + l)"+5 (a;j/)~^<''+^^"^ [(i - x)(l - 
X [xy(l — a;)(l — y)] ^ cos cos Ay 



A^ [n{a + 1) + to + ]^]e^ + [n(6 + 1) - to + ^]7:r - (n + i)/?:, + J 



(6.19) := [n(a + 1) + + + [n{b + 1) - so + i]7a ^ {n + ^)py + J. 

We are justified in ignoring the terms 0{l/p) in (16. 6p due to it being a uniform bound. Combining 
([eTf]) and ([618]) shows 

(6.20) c„Fi''"+*°'''"-*°)(a;)P^"'+'*«'''"-"«)(y) 



n'*°~'°(a + 1)^^(6+ 1) 



1 



cos (Aj:) cos (By) 



Next we make use of the trigonometric identity cos (A) cos {B) — ^ (cos {A + B)+cos {A — B)). 
The term cos {A + B) oscillates with frequency proportional to p so contributes at a lower order 
to the summation than the term cos{A — B) which has argument of order unity. Now from (|6.14p 
and (Pl^ 

(6.21) cos(A, - Ay) cos {w{X - Y){e,{a + 1) + 71(6 + 1) - pi) - {sq - to){eo - 70)) 
which using (|6.15p can be written 



(6.22) 



cos(Aa; — Ay) ^ Re exp 



twiX - y)^/^ 



[ AuX,{l~Xs) 



„'i(so-io)(eo-7o) 
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Further, we have from (|6.15p that 
(6.23) 



gi(so-to)(So-7o) 



Substitutmg in (j^?^ . multiplying by 1/2 and substituting in (jOO)) gives (j6.16p . The 

summand for < w < 1 contributes to a lower order because with A > the Jacobi polynomials 
are exponentially small [20l [16] . □ 



This result in turn allows the leading large p behaviour of the sum in (j6.10p and thus the limit 
of the scaled correlation kernel (|6.2p to be computed. 



Proposition 6.4. Let K{so,Y; to, X) be as specified by i6.^) . and let S be such that I < S < fc+1 
so that r{pS) — p. Introduce too the notation 



2+k 1+k 
(6.24) 1/ := ' 



k \IS{2 + k-S)' 
We have 

(6.25) K{so, Y; t^, X) = jr (sq, to, X) 
where 

(6.26) A = e-, {2 ^ kfkSi2 - k - S) 



and 




4 + 8fc + P(l + S) + /c2(5 + 25 - 52) 
Re (e"*(-^-^)(l + itz/)""-*") dt, so > to. 



(6.27) K*iso,Y;to,X). 
Proof. It follows from Proposition [63] that 

„_i (S-l)p (l + fc-S)p 

(6.28) ^ c„P^'^"+*°'^"-*«)(a;)P^"'+"° )(y) - -^—^ ( ^ ) ( f ) 2: 

where, with wo such that — Aq = 0, 

dw 



sq — U-) 


j dt, 


So 


< to- 




(S-l)p 








2 


(r 













ttV-Ao 



iw{X - y)V^ \ / 2aw(l - Xs) + 26wXs + iwV^A^ 
''''P iv 4ziXs(l - Xs) ) \ AXs{l-Xs) 



Recalling the dependence of a and 6 on ui from ()6.12p . and making the substitution 



47rXs(l - Xs)us 
gives 



(6.29) 1 = J^^Re 



exp{TTit{X - Y)) 



(2 + fc)2(fc2g ^ 2iGt + fc(25 - 52 + iGt)) \ 
4 + 8fc + /c3(l + 5)+A;2(5 + 25-52) j 



where G = ^5(2 + fc - 5)(1 + /c). Substituting ([05)1 in (jg?^ and recalling (jgTUD we obtain, 
after some further minor manipulations, the result (|6.25p in the case sq > ^o- As noted below 



24 



BENJAMIN J. FLEMING, PETER J. FORRESTER, AND ERIC NORDENSTAM 



(I6.12p . in the case sq < to, instead of (|6.10p we have 

oo 

K{pS + so,y;pS + to,x) ^ -Y,Cn{x,y)Pt^+'"''"-'"Hx)Pi''''+'"''>"-'»^ 

n—p 

Thus, up to the minus sign, the asymptotic form is again given by (|6.28p . but with the terminal 
of integration in X now from 1 to cx). This gives ()6.25p for sq > to- 



Corollary 6.5. With the scaled correlation function p specifed by !16.4\ ), and K* specified by 
we have 

(6.30) p{ti,Xi; . . . ; X„) = det[K*{X„ tf, X,, t,)],j=i,...,„. 

Proof. In the region 1 < »S' < fc+l, this is an immediate consequence of Proposition [6^ (note that 
the prefactors ^^^^S*o^*o (j6,25p cancel out of the determinant). The correlation function 
in the region l + fc<S'<2 + fc follows from the form in the region < 5 < 1 upon making the 
mappings 

{X,,t,)^ i~X„~t,), S^2 + k-S 

which from (|6.27p is indeed a symmetry of (|6.4p . To show that (|6.4p is valid for < 5* < 1 
requires deriving the analogue ()6.4p in this case. Recalling I5.55P and ()5.59p . and with sq > to, 
we must then obtain the large p form of 

(6.31) K{pS + so, y;pS + to, x) = £ y)PLtl;^f''''''\x)PLtl7n^^^^^^ iv)- 

n=0 

We have done this, using the same general strategy as for 1 + fc < S* < 2 + fc, obtaining a result 
consistent with (16.41). □ 



In the Introduction it was remarked that the bead process was introduced by Boutillier [H] 
as a continuum limit of a dimer model on the honeycomb lattice. The corresponding scaled 
correlation was calculated to be of the form (|6.30p but with K* replaced by , where 



(6.32) J^iso,Y;to,X) 



]- J'^ (e'*(^-^)(7 + *iv/r^)^«-*") dt, so > to, 

J (e'*(^-^H7 + ^tVl^y°~'°) dt, So < to- 

JR [-1,11 ^ ^ 



The parameter 7, I7I < 1, represents an anisotropy in the underlying a6c-hexagon, with 7 = 
corresponding to the symmetrical case. 

We observe that the factor (7 + ity^l — j'^yo-ta ^j^g integrands of (|6.32p can be replaced 
by (1 + ity^l — 72/^)so-*o without changing the value of the determinant. Changing scale 
J^{so,Y;to, X) 1-^ nJ^{so,TTY;to,TTX), and comparing the resulting form of (|6.32p with (|6.27p 
shows the two results to be the same, upon the identification 



(6.33) yW_2 + M 1 + 



7 k \]S{2 + k-S)' 

although the quantity on the RHS is always positive. This latter feature is a consequence of the 
calculations relating to our finitized bead process being carried out under the assumption that 
q > p. From the symmetry of the hexagon, the case q < p is obtained by simply replacing X, Y 
in (|6.32p by —X, —Y, or equivalently replacing i' in (|6.27p by —ly. This then allows us to extend 
(|6.33p to the region —1 < 7 < 0, by replacing the 7 in the denominator by |7|. 
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Appendix 

Chen and Ismail [2^ use the asymptotic method of Darboux (see e.g. [24 ) applied to the 
generating function 

oo 

(A.l) = (1 + + V)''+'[1 - < - 6ry - (1 + a + b)^7j]-^ f{t), 

where ^ and rj depend on x and t according to 

(A.2) ^=k^ + 1)^(1 + O'+ni + and = 

2 X + 1 

The basic idea is to identify and analyze the neighbourhood of the i-singularities of f{t), to 
replace f{t) in (jA.ip by its leading asymptotic form g{t) in the neighbourhood of the singularities 
(referred to as the comparison function), and finally to expand the latter about the origin to 
equate coefficients of and so read off the asymptotic form of P^^°-^'^^^^{x). 
It is shown in [20] that the t-singularities of f{t) occur at 

/A Q^ ^ b{x-l) + a{x + l)±l^/^ i_ 

where A is given by (|6.5p and 

2(1 + a + oj(l — a; + 1 

The comparison function is computed as 

(A.5) g{t) = B+{t+ - t)-'/' + B^ih - t)-'/' 

where 

(A.6) B±= lim(t±-t)i/2/(t). 

t — n± 



Since we know from [20] that 
dt 

we calculate from (|A.6p that 



dH 



= : 2A±, 



(l+x)2(l+e±)2+a(l+^^)2+b 



„ ^ (l+e±)°+^(l + r?+)^+H^ + l)v/^ 

(A. 7) B±^T . , yr- 

= eT"'/4(-A)-i/4(l + e±)"""/'(l + r/±)^-''/2 

The first of the formulas in (|A.7[) is not reported in [5D], while the second is their (2.14), (2.15) 
but with our e^^^^^{~A)^^^^ replaced by —i{A)~^^^ and i{—^/A)~^/^ respectively. 
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The coefficient oft" in (|A.5p . and thus the leading asymptotic form of p^+a^./S+^n^-j;) according 
to the method of Darboux, is equal to 



V» I 2 

n 



(A.8) (-1) 
To simplify this, we note from (|6.5p and (|A.4[) that 

1 + C± 



[B+t" +B^t_'^ ^] 



{1- x){a + b+l) 



l + V± = 
2e± 



2(&+l) V^\±r 



(1 +x)(a + 6+ 1) 



+ l ^{a + b+l){l-x^) 
These substituted in (|A.3p and (|A.7p give 



B.f 



1 



2 cos /i(6', 7, p) 
2(6+1) 



2(a + l) 



(l-a;)(l + a + 6) 
1 t(f'+i)+l+3 



Ka+l) + f + i 



(l + a;)(l + a + 6) 



(l-a;^)(a + &+l) 



where 



h{e, 7, p) = [n(a + l) + a + ^]9+ [n{b + 1) + /? + i]i7 - {n + + J 



This, together with the expansion 



substituted in (fO)) gives (p^ . 



(-1) 



n ( 2 

71 



n 2 
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